Ostwald Ripening in Two Dimensions: Treatment with Pairwise Interactions. 
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We present a systematic extension of mean field models for Ostwald ripening in two dimensions. 
We derived, using a mean field type approximation, an analytic expression for pairwise interac- 
tions between the minority phase droplets. These interactions appear in dynamic equations for 
the droplets' radii and positions, which are solved numerically for systems of tens of thousands of 
droplets, yielding very good agreement with recent experiments. 

PACS numbers; 64:60My 64.60Cn, 64.75+g 



Ostwald ripeningp is a coarsening process, observed 
during the late stage of the evolution of a two-phase 
system (say, solid in liquid), in the course of which 
the droplets of the minority phase exchange material by 
means of diffusion. This process leads towards a scaling 
state in which the characteristic length scale grows with 
time according to the Lifshitz-Slyozov lawo R{t) ~ t-^^^. 
When rescaled by R{t), all statistical characteristics of 
the system (such as droplet size distribution, spatial cor- 
relations etc) are time-independent. 

These characteristics were studied in a number of 
detailed |-|iiumerical simulations. The Cahn-Hilliard 
equatiorJja, solved in the entire system, provides the 
most basic model of Ostwald ripening. If the solid 
droplets are assumed to be uniform and only their bound- 
aries are retained to describe the system, the problem 
becomes one of solving the diffusion equation for the 
concentration field c(r) between droplets, with Gibbs- 
Thomson boundary conditionsO: 



c\droplet{t) = Ceq{R{t)) = Coo + 



R{t)- 



(1) 



The boundaries move in response to the incident flux; 

J = - Vclfl (2) 



dR 



At long times a quasi-static approximation can be used 
and the diffusion problem turns into Poisson's equation 
in 2-d, 



V'c = 0. 



(3) 



This problem can be reduced to an implicit system of or- 
dinary differential equationaJa. Since these calculations 
take into account all the complicated interactions be- 
tween the droplets, mediated by the diffusion field, they 
do not elucidate the relative importance of different as- 
pects of these interactions. The number of droplets and 
the length of the simulations are also limited. On the 
other end crfiUiei complexity scale are analytical mean field 
treatmentsaa^Ei] that neglect all spatial effects; these, 
however, are too simple to account for correlations in 



a system of non-vanishing volume fraction. Previous at- 
tempts to "interpolate" between these extremes, and|prp=. 
vide an analytical treatment of the spatial effectalj'E£l 
contain uncontrolled approximations and lead to very 
complicated expressions. 

Using a mean-field type approximation, we calculated 
analytically a pairwise interactions between the droplets. 
These appear in closed-form dynamic equations for many 
interacting droplets which are integrated numerically. 
Our model gives rise to a very efficient numerical al- 
gorithm; the evolution of tens of thousands of droplets 
can be followed. Our work was motivated by and our 
results are compared with a recent experiment on a two- 
dimensional fika of liquid and crystalline succinonitrile 
in coexistencell3. This experiment shows that even at 
if = 0.4 the droplets are almost circular; therefore we 
characterize the system by the set of the droplets' radii 
Ri and the positions fi of their centers. 

We write the solution of Eq. (^) in terms of " charges" 
and "dipoles" placed at fi 



{p^ ■ (f- fj)) 



(4) 



(i?o is an arbitrary length). The boundary conditions (|^) 
on the surface of each droplet give rise to linear equations 
which determine the charges qi and the dipoles pi. When 
the resulting solution is used in (^) to derive the flux at 
droplet i, we find that it's normal component has two 
terms. One is isotropic, due to the charges; it affects the 
droplet 's area and therefore determines the dynamics of 
the radii. The other is an anisotropic contribution, due to 
the dipale.s, which gives rise to a shift of the centers of the 
dropletaUJo. The charges and dipoles can be eliminated, 
yicldingla dynamic equations for the radii and positions 
of the droplets: 



R^ 



qi_ 
Ri 



'dt 



^Lrj(l + i?,-/i?,) = ^Lrj, 



R? 



2^RjR.j 



\X^, 



(5) 



(6) 
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where Xi^ = \fj — fi\ and Rc = a/coc is a capillary 
lenffitijp-| The second equality in Eq.(|5|) is due to a sum 



rulel 



are the elements of the inverse of the ma- 



trix L, defined aa 



Li'j — RiRj 



log(i?,/i?o) 
log(X,,,/i?o) 



if i = j 
otherwise 



(7) 



The value of the parameter i?o is axhitrary; it does not 
affect the dynamics of the systemEJ. Once L has been 
inverted, we can evaluate the radii at the next time step. 
One must invert the N x N matrix Lij at each time step 
(which takes 0{N^) operations); since we need 0(N) time 
steps, each rua would cost 0{N'^) computations. 

Beenakkeilij solved an analogous problem in 3-d by 
truncating,,the matrix Lij. Motivated by the effect of 
screeningLIx3, he took into account only interactions be- 
tween those droplets whose separation does not exceed a 
threshold. Akaiwa and MeironEl used the analogous trun- 
cation procedure in 2-d. However, formal truncation of 
the matrix seems problematic. Since the matrix elements 
grow with Xij, the elements of the inverse matrix 
as functions of the cutoff should contain large fast oscil- 
lating comunnents. Apparently, the success of Akaiwa 
and Meironia is due to the fact that these oscillations are 
effectively averaged out during the run. 

Rather than discarding large matrix elements and pro- 
ceeding to invert the truncated matrix numerically, we 
calculate analytically the elements of L~^^ using a mean 
field approximation. First, represent L as the sum 
of its diagonal and off-diagonal parts, L = Lq — Li, 



where (Lo)jj = Si^Rf log{R^/ Rq) and (Li) 



Sij)RiRj log{Xij / Rq) and introduce a new matrix T de- 
fined by 



^0 ^^-^0 ^ 



This is indeed the inverse of L if T satisfies the equation 



T ^ Li+LiL^^T 



(8) 



For convenience we represent Tij — RiRj(f)i,j; Eq.(|8|) now 
becomes a condition on the matrix (jj, which we write 
separately for the diagonal and off-diagonal elements: 



log(Xfc.j/i?o) 



[ log{R,/Ro) 



(9) 



E t!^S^'^.-^^ + '^^'^ = log(X,,,./i?o) (10) 



with 



log(i?,/i?o) 



7fc = 1 



4>k,k 



logiRk/Ro)- ^^^^ 

The set of equations (|^^l|) is exact but solving it is as 
difficult as inverting L. The advantage of this formalism 



is that it serves as a convenient starting point to gen- 
erate approximate expressions for - in particular, a 
manageable mean-field approximation. 

We give now an outline of the main ideas and formal 



steps of the derivation of our approximation to 



^2 J I 



de- 



tails will be presented elsewhereo. First of all we assume 
that the i,j matrix: element depends only on the corre- 
sponding distancellj, = 4>{Xij), and that (j){X) varies 
only on scales ^ such that the number of droplets within 
an area is large. If this holds we can replace in Eqs.(|lO|) 
and (||) the factor 1/ log(i?j7i?o) by its mean value and 
the sums by integrals. Within this mean-field approxi- 
mation Eq.(^ imphes (/)[™'^2_^=const. = </>o; this leads to 
a self-consistent expressions for (j){Xij) only if we set 
7fe — const. = 7o which, in turn, means (see Eq.(|ll|)) that 
(j)Q — and 7o = I. Thus within our mean-field treatment 
of 0, Eq.([lO|) becomes 

2^ /log(^)<^(%fe)d'r,-0(X,,) = log^ (12) 
where we introduced a constant length, defined by 



= 27m 



1 



\ogiRo/R) 



2iTn 



log(i?o/ (i?)) 



(13) 



Here n ~ ip/Tr{R'^) is the number of droplets per unit 
area; the angular brackets denote averaging over the 
droplet size distribution. Operating with on Eq.(|T^) 
yields the differential equation 



r'0 + v'0--27r5(f). 



(14) 



which is solved by the zeroth order modified Bessel func- 
tion of the second kind (also called MacDonald function) , 
4>{X) = Ko{X/(). This determines the off — diagonal 
elements of (jjij up to the length C which depends on the 
arbitrary parameter Rq. Since the value of this param- 
eter does not affect the exact solution of the problem, 
we can tune it to improve our approximation. Note that 
when our (t>{X) is used in the integral form of Eq.(^ we 
get 

^ J \og{Xk,j/Ro)KoiX,^k/Od^r, = cj^k^k - <^o - 

(15) 

After some lengthy algebra, presented elsewhere0, we 
can show that this condition is satisfied when Rq satisfies 
the following approximate expression: 



log(i?o/i?) 



Ko(R/o) 



(16) 



By substituting Eq. ([iq ) in Eq. ( |l3| ) we eliminate Rq from 
the problem and obtain an equation for (: 



2Tm 



1 



\Ko{R/0/' 



(17) 
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which is almost identical to Marqusee's expression for the 
screening lengthEJ. The solution for (j){X) determines 
the matrix T and we get an expression for the matrix 
elements which, when used in Eq.(||) together with 
the relatiorO 

finally yields our central result: 

^ Koix,./C) (I I] .3. 

KoiR,/0 R^J^ ^ ' 

Since ^^0(2;) ^ exp(— a:) for large x, C, indeed has the 
meaning of a screening length. Note that the total area 
of the droplets is conserved by Eq.(|l8|). 

The approximations made, of replacing the sums by in- 
tegrals and l/log(i?j/i?o) by its average value 1/{2tt('^), 
are valid if the function (fi{X) is smooth and the number 
of droplets in the screening zone is large: iV^ — nC,"^ ^ 1. 
For if = 0.13 (as used in the experiment) the approach 
outlined above gave a much too small value for A 
somewhat larger value was obtainedEj by taking into ac- 
count the fact that each of the droplets is surrounded by 
a depletion zone, from which all neighbors are excluded. 
Including de plet ion zones retains an equation of the same 
form as Eq.([18|) for the evolution of the radii, but with 
the mean-field ( replaced by > C- For V = 0-13 
it becomes (sc — 2.73R « 0.56X where X is the typi- 
cal distance between neighbors. Although (sc is still too 
small to provide a formal support to our approach, we 
tested it by integrating (numerically) Eqs.([l8|) and (g). 

The detailed descpiption of our numerical algorithm is 
presented elsewheretil. An advantage of our method is 
that it conserves the total area of the droplets at each 
time step. Computational efficiency is very high: an 
entire run, starting with N droplets initially and run- 
ning till most disappear (i.e. reach Ri « 0) takes only 
O(A^logA^) operations. This is due to the fact that we 
are able to keep the time step reasonably large, eliminat- 
ing a large number of droplete-at-each step. At the same 
time, unlike previous studiesEl'll3'E3, we ensure exact con- 
servation of the total area of the droplets at each time 
step. This makes our approach useful for an extensive 
study of the Ostwald problem. 

We used initial systems of 50000 droplets with toroidal 
boundary conditions. First evaluate Ri using Eq.dl^) 
with C — 2.73R; next use Eq.(||) to calculate the motion 
of the droplets' centers. 

Regarding this mption, one can use Eq.(|6|) to esti- 
mate its importancetllEJ by comparing the characteris- 
tic times for significant shift versus a shrinking droplet 's 

v^ 

lifetime; Tshift/TiifG '~ ~ T^/f- Hence for small frac- 
tions (and even for f = 0.13) the shift of the droplets is 
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FIG. 1. The size distribution of the droplets; R is nor- 
malized by R. Data were taken when the number of droplets 
(initially 50 000) was iV = 3000 {largest circles), N = 2000 
{medium circles and N — 1000 {small circles. Each plot 
presents the average over 8 runs; the solid line is a guide 
to the eye. The fact that at these three times the same dis- 
tribution was obtained indicates that we reached the scaling 
state. The experimental result of Krichevsky and Stavans is 
also shown {diamonds). 

adiabatically slower than the droplets' growth and one 
is tempted to neglect the motion of the centers. Indeed 
we tried this, but the resulting solution of the dynam- 
ics of the radii alone led to severe inconsistencies (man- 
ifest in significant overlap of neighboring droplets, i.e. 
Ri + Rj > Xi,j). This convinced us that the correla- 
tions induced by the proper movement of the droplets 
are essential for a physically sensible solution. 

The sum on the r.h.s. of Eq.(|^) appears to have con- 
vergence problems; assuming that the charges qj are un- 
correlated random variables, with zero mean (as required 
by the total area conservation) one can easily get 

where Lg is the size of the system. The charges, however, 
are correlated, (the effect of screening) and this causes fast 
convergencdiJ. We performed in some cases the entire 
sum and then repeated the simulations, summing only 
over droplets j with \Xij\ < b, using b = 2.3X and b = 
3.25X (the latter value corresponds to about 10 terms 
in the sum) ; no significant difference between these three 
simulations was found. 

For 50,000 initial droplets at = 0.13 the sys- 
tem approached the scaling state at iV « 3000; this 
is when we started our measurements. In order to 
reduce fluctuations we averaged all the data over 8 
runs. Such a run took about 105 minutes on a sin- 
gle processor of an IIP-9000 (series K200) worksta- 
tion. First we present the distribution of the droplets' 
radii in the scaling state (see Fig. J). The agree- 
ment with the experimental pointsEj is reasonable. 
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FIG. 2. Weighted time average of the correlation func- 
tions of droplets' positions, G{r), taken at different times in 
the scaling state, each averaged over 8 runs {circles); the ex- 
perimental results of Krichevsky and Stavans (squares). The 
solid line is a guide to the eye. 

A much more demanding comparison is that of various 
correlation functions Fig. 2 presents the correlation func- 
tion G{r) - the probability to find a droplet's center at 
a distance r from a given droplet's centeiu. as obtained 
by our numerical solution and experimentlij. We recover 
the positions of the initial rise of G{r) as well as the ex- 
istence and position of its maximum. The value of our 
curve at the maximum is smaller than the experimental 
one; our G{r) is rather close to the result of Masbaum, 
obtainedJpy solving the Cahn-Hilliard equation (see Fig. 
la in ref.Q), that also exhibits a small peak. 

We measured also the— correlation functions of the 
charges defined first in ref.E£l, which contain more detailed 
information about the system. For a charge qi calculate 
Q+(r), the total amount of similar charge as qi within an 
annulus [r, r + dr] around and define the function 

9+{r) = {qiQ+{r))- 

Similarly we define g-{r) in terms of the opposite charges. 
These two functions, as obtained by our simulations, are 
presented in Fig. 3 together with the corresponding ex- 
perimental data of Krichevsky and Stavans. The agree- 
ment between theory and experiment is quite impressive. 

In summary, we introduced and tested a model for 
the Ostwald ripening process in two dimensions. Our 
model approximates the exchange of material between 
the droplets, as determined by the complicated diffu- 
sive interaction, by simple pairwise couplings. An ap- 
proximate mean field type approach is used to evaluate 
these pairwise interactions. The resulting dynamic equa- 
tions are then solved numerically. Our model gives good 
agreement with experiment for a fairly large value of area 
fraction. We found that the shift of the droplets plays 
an important role at ip — 0.13. Our method leads to a 
very efficient numerical algorithm that can be useful for 
future studies of this problem. This research was sup- 



FIG. 3. The correlation functions for the same and op- 
posite charges, as obtained by our simulations and by experi- 
ments. Our data present averages over 8 runs. The solid line 
is a guide to the eye. 
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